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Uses of Thermographic Simulations 

Description of Flash Thermographic Measurement 

Thermal Simulation of Heat Diffusion in Unidirectional and Quasi- 
isotropic composites 

Simulating delaminations in composites with quadrupole method 



NASA Applications of Thermographic 
Simulations 



• Optimization of thermographic techniques 

• Determination of limitation of experimental and analysis techniques 

• Generation of set of eigenvector of principle component analysis of 
thermal response 

• Inversion of thermal responses to information on flaw sizes and 
locations 

• Training neural networks for rapid reduction of thermal data 




Flash Thermography Measurements 


Specimen 


Flash Heat Source 


IR Filter 


Light Shield 

Infrared Camera 


Data Collection : 

TWI Echotherm System 
FUR SC6000® infrared imager 
640 x 512 element array 
12"x9" = ~0.75 ft 2 


3-5 pm wavelength range 
120 fps acquisition speed 







Simulation of Composites 


Finite Difference 

• Implicit and explicit 

• Initial flux at surface or initial temperature 

• "Finite Difference Methods in Heat Transfer" M. N. Ozisik 

Finite Element 

• Simpler for complex geometries 

• Initial flux at surface or initial temperature 

• Commercial packages are available 

Thermal Quadrupoles (examined in this paper) 

• Modification of method presented in "Thermal Quadrupoles, Solving the Heat 
Equation through Integral Transforms", Maillet, Andre, Bastsale, Degiovanni and 
Moyne 



Thermal Response to an Impulse Flux Input at 
Single Point (S(x) 8(y) 8(z) 8(t)) Unidirectional 
Composite ( a x = 10 a ) - 1/60 sec after pulse 



Front surface temperature Back surface temperature 

Front and back surface temperature are approximately the same 
Fiber direction evident in thermal response 

Single Ply 0.01 cm thick 





Thermal Response to 8(x) 8(y) S(z) 8(t) Quasi Isotropic 
Composite with Two Plies, Top ply ( a x = 10 a A Second ply( 
10 a x = a ) - 1/60 sec after pulse 


Front surface temperature 

■ Two Plies, both 

0.01 cm thick 


Back surface temperature 



• Initial front surface temperature indicates upper ply fiber direction and back surface 
temperature has appearance of isotropic in-plane heat flow. 

• Back surface response is more indicative of flaw response 
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Gaussian Fit of Temperature Profiles in 
Simulated Responses as Function of Time 


Front surface temperature profiles Back surface temperature profiles 



For quasi-isotropic layup the in-plane diffusivity can be considered to be approximately isotropic 


Composite Flat Bottom Hole Specimen 



Approximate Depths of Flat Bottom Holes 0.05, 0.1, 0.15 cm 
Approximate Diameters 1.27, 0.63, 0.32 cm 
Quasi-isotropic ply layup 



Temperture (au) 



Temperature Profiles in X and Y directions over Holes in 
Composite Specimen With FEM Simulation Assuming 
Isotropic In-plane Heat Conduction 




0.32 cm diameter hole 0.05 cm below surface 


0.64 cm diameter hole 0.1 cm below surface 


Thermal Quadrupoles for Single Layer 
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• T 0 - Laplace transform of front 
surface temperature 

• F 0 -Laplace transform of front 
surface flux 

• T r Laplace transform of back 
surface temperature 

• Fj -Laplace transform of back 
surface flux 

•a - Thermal Diffusivity 

• K - Thermal Conductivity 

• d - layer thickness 


Two of T 0 ,F 0 ,T 1 ,F 1 are defined, inverse transform solved analytically 





One Dimensional Thermal Quadrupoles for 
Two Layers 
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• T () - Laplace transform of front 
surface temperature 

• F 0 -Laplace transform of front 
surface flux 

• T 2 - Laplace transform of back 
surface temperature 

• F 2 -Laplace transform of back 
surface flux 

• a.|,a 2 - Thermal Diffusivities 

• K 1? K 2 - Thermal Conductivities 

• d 1? d 2 - layer thicknesses 


Two of T 0 ,F 0 ,T 2 ,F 2 are defined, inverse transform solved numerically 





Laplace Transform Thermal Response in 
Plate 
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K z , diffusivities in x, y, z directions and thermal conductivity in z 
Cosine Fourier coefficients for surface temperature and flux 



Two Dimensional Thermal Quadruples for Two 
Layers with Zero Flux at Back Surface and No Spatia 
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Variation in Front Surface Flux(F 0 ) 

T 0 (x,s),F 0 

• T 0 , T la , T lb , T 2 - Laplace transform of 
temperature at front surface, above 
and oelow interface and back 
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• Fq, Fj, /s- Lap lace transform of 
front surface interface and back 
surface flux 

• a b a 2 - Thermal Diffusivities 

• K b K 2 - Thermal Conductivities 

• dj,d 2 - layer thicknesses 

• m- refers to Fourier cosine series 
coefficient 


T 0 , T la , T lb , T 2 can be defined in terms of F 0 and Fj 




Solving for Temperature at Front Surface 


Instead of solving for Fourier cosine series coefficients, discretize the 
temperature and fluxes at surfaces and interface 

Tob(xn), T la (x n ), T lh (x n ) , T 2 (x n ), and F 1 (x n ) 

• Represent the all temperatures and fluxes as vectors, then relate temperatures to 
fluxes with matrix equations 
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Delamination Response (au) 



Estimating the Flux at the Interface as Zero Flux over 
Delamination and 1 Dimensional Flux Value in Region with 
No Delamination 




Position (cm) Position (cm) 


Comparison of FEM solution to approximate quadrupole 
solution 



Setting Temperature Equal Below and 
Above Interface 
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Resulting in the matrix equation 
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Matrix equation solved numerical 
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Delamination Response (au) 



Temperatures and Fluxes Calculated for 
Interfaces and Front Surface 



Quadrupole solution calculated in 0.3 sec, FEM solution 
in 30 sec 



Summary 



• Simulations assuming in-plane thermal conductivity is isotropic is a 
good approximation of quasi-isotropic composite layups 


Quadrupole method is a computationally efficient technique for 
simulating the thermal response of delaminations in composites 


